Probing the existence and dynamics of Majorana fermion via transport through a quantum dot 
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We consider an experimentally feasible setup to demonstrate the existence and coherent dynamics of Ma- 
jorana fermion. The transport setup consists of a quantum dot and a tunnel-coupled semiconductor nanowire 
which is anticipated to generate Majorana excitations under some conditions. For transport under finite bias 
voltage, we find that a subtraction of the source and drain currents can expose the essential feature of the Ma- 
jorana fermion, including the zero-energy nature by gate-voltage modulating the dot level. Moreover, coherent 
oscillating dynamics of the Majorana fermion between the nanowire and the quantum dot is reflected in the 
shot noise, via a spectral dip together with a pronounced zero-frequency noise enhancement effect. Important 
parameters, e.g., for the Majorana's mutual interaction and its coupling to the quantum dot, can be extracted out 
in experiment using the derived analytic results. 
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The Majorana fermions, proposed in 1937 by Majorana 
[|lj], are exotic particles since each Majorana fermion is 
its own antiparticle [|[. Recently, the search for Majo- 
rana fermions in solid states, as emerged quasiparticles 
(elementary excitations), is attracting a great deal of at- 
tention In addition to being of interest for some 
fundamental physics such as the non-Abelian statistics, 
the Majorana fermions emerged in solid states are also 
a promising candidate for fault-tolerant topological quan- 
tum computation (TQC) [JTcj— [T3|] . In solid states, the pre- 
dicted Majorana bound states (MBSs) are coherent super- 
positions of electron and hole excitations of zero energy, 
which can appear for instance in the 5/2 fractional quan- 
tum Hall system[|l4|| and the p-wave superconductor and 
superfluidnla]. In practice, an effective p-wave supercon- 
ductor can be realized by a semiconductor nanowire with 
Rashba spin-orbit interaction and Zeeman splitting, and in 
proximity to an s-wave superconductor . 

Of crucial importance is then a full experimental obser- 
vation of the Majorana fermion in solid states, from multi- 
angles in a variety of contexts to demonstrate its existence 
and coherent dynamics. Some existing proposals include 
analyzing the tunneling spectroscopy which may reveal 
characteristic zero-bias conductance peak [|l6| and pecu- 
liar noise behavior s jl7[ |l8| ], verifying the nonlocality na- 
ture of the MBSs [|19|, [20|], and observing the 4-tt periodic 
Majorana-Josephson currents [[I). We notice that 

a breakthrough was achieved in a most latest experiment 
[ f22| l where, using the tunneling spectroscopy, a mid-gap 
Majorana state was detected, based on the observation that 
a zero-bias-peak stickily appears in the differential con- 
ductance through an InSb nanowire, while the nanowire is 
in contact with a superconducting electrode (NbTiN). 

In this work we consider a setup based on trans- 
port through a quantum dot (QD) to probe the existence 
and coherent dynamics of Majorana fermion, by weakly 



tunnel-coupling the QD to a nanowire as realized in the 
above-mentioned experiment [p2[]. We find that a sub- 
traction of the source and drain currents can expose the 
essential feature of the Majorana fermion, and the zero- 
energy nature can be justified by gate-voltage modulating 
the QD level. Moreover, oscillating dynamics of the Ma- 
jorana fermion between the nanowire and the quantum dot 
can be identified from the shot noise, via a spectral dip to- 
gether with a pronounced zero-frequency noise enhance- 
ment effect. 
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FIG. 1 : (Color online) Schematic setup for detecting Majorana 
fermion by transport through a quantum dot, while the dot is 
coupled to a semiconductor nanowire which is in contact with 
an s-wave superconductor. Under appropriate conditions such 
as with large Zeeman splitting and strong spin-orbit interaction 
(SOI), a pair of Majorana bound states is anticipated to emerge 
at the ends of the nanowire. 

Model. — Figure 1 shows schematically the setup that 
will be studied in this work, where the transport through 
the quantum dot (QD) is employed to detect the Majo- 
rana fermion emerged in the tunnel-coupled nanowire. 
Combining a strong Rashba spin-orbit interaction and the 
Zeeman splitting, it was shown that the proximity-effect- 
induced s-wave superconductivity in the nanowire can 
support electron-hole quasiparticle excitations of Majo- 
rana bound states (MBSs) at the ends of the nanowire [g- 
[7]]. Since the Zeeman splitting should be large enough 
in order to drive the wire into a topological supercon- 
ducting phase, we can assume it much larger than the 
transport bias voltage and the tunnel coupling amplitude 
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and rates. In this case, we can model the QD by a sin- 
gle resonant level and treat the electron as spinless par- 
ticle. Accordingly, the whole system of Fig. 1 can be 
modeled by H = H Lc&ds + H t 



sys 



Leads 



J2 a =L,RJ2k £ akcl k c ak describes the normal metallic 

leads, H T = J2 u =l,r J2k t akdc^ ak + H.c. is for the tun- 
neling between the leads and the dot, and the low-energy 
effective Hamiltonian for the central system reads 



e_ D rf t (i+ -e M 7i72 + (Ad- A*d t )7i. (1) 



Here c^ fc (c a fc) and ct>(d) are, respectively, the electron 
creation (annihilation) operators of the leads and the dot, 
with the corresponding energies e ak and ej> Particularly, 
in Eq. the second term describes the paired MBSs gen- 
erated locally at the ends of the nanowire but with a mu- 
tual coupling 6m ~ e~ £ / ? , where £ is the wire length and 
£ is the superconducting coherent length. The last term 
in Eq. ([[]) describes the tunnel coupling between the dot 
and the left MBS. For spinless dot level, we can choose a 
real constant for A, while in general A has a phase factor 
associated with the spin direction. 

To solve the transport problem associated with Eq. ([[]), 
it is convenient to switch from the Majorana representa- 
tion to the regular fermion one, through the exact trans- 
formation 71 = f< + f and 72 = i(f' — /). / is the 
regular fermion operator, satisfying the anti-commutative 
relation {/, f'} = 1. Accordingly, we rewrite _ff sys as 
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+ A(d-dt)(/t+/). 

(2) 

For latter convenience of discussion, we rearrange the tun- 
nel coupling term in Eq. (^) as Hi = (Xf^d + Xipdf) + 
H.c. where Ai = A or corresponds to coupling to a MBS 
or a regular fermion bound state. In the transformed rep- 
resentation, the basis states of the central system are given 
by \ndUf), with rid and rif being or 1 so that we have 
four basis states {|00), |10), |01), |11}}. 

Rather than the linear response[^6[, we consider trans- 
port through the quantum dot under finite (large) bias 
(in this regime the temperature can be approximated to 
zero in calculating the transport currents). This case cor- 
responds to the single dot level deeply buried into the 
voltage window, which allows us to employ the standard 
Born-Markov master equation to solve the problem. For 
the transport setup, the two leads are regarded as a gener- 
alized environment, and the embedded system in between 
the leads as the system of interest which is described by 
the reduced density matrix p(t). For the latter use of cal- 
culating the noise spectrum, we present below the particle- 
number resolved version of the master equation for p( n > (t) 
as pi 



(dUp (n) + p {n) dU - 2V"~ 1) d t ) , (3) 



with "n" the electron number transferred through the cen- 
tral system and satisfying Y^=o P^ W = P(t)- Here 
we introduced the Liouvillian as Cp = [H sys , p], and the 
tunneling rate of the dot to the lead as T a = 2Trg a \t a \ 2 , 
with g a the density of states of the lead a (L or R). Cor- 
responding to Eq. ([|), the unconditional Lindblad master 
equation reads p = — iCp + T iX)\d^]p + Ti{D[d]p, where 
V[A]p = ApA^ - ±{A^A,p}. 

Current. — First we consider the transport current, 
which can be evaluated in a rather simple way as follows. 
Instead of solving Eq. (^|), we evaluate the unconditional 
Lindblad master equation for p(t), and calculate the aver- 
age occupation of the quantum dot, fid = Tr [c^c^i)] . 
Then, the left- and right-lead currents are obtained as 

hit) =eT L (l-n d ), I R (t) = eT R n d . (4) 

For steady-state currents, analytic results read 
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Here we introduced the total rate V = Tl + Tr, and the 
asymmetric factor y = Tr/Tl. Also, in the above results 
we restore Ai = A as it should be for coupling to Majorana 
fermion. 

We notice that the first term, eTiXn/T, is just the 
steady-state current through a single resonant level with- 
out any side coupling, or coupling to a regular fermion 
bound state. In the latter case, the side coupling does not 
affect the stationary state, but only affects the transient 
behaviors (e.g., the noise properties). The second term, 
in the ^(ii)' originates from the Majorana nature, partic- 
ularly its hole excitation (antiparticle component). How- 
ever, of interest is that this term appears as a consequence 
of the interplay of the electron and hole excitations in the 
Majorana fermion, but not their separate individual role. 

Moreover, an unusual feature appears that, even in 
steady state, the source and drain currents are unequal to 
one another, with a difference of 



A/ = I R -P L = 



4e(T R - T L )X 2 



r 2 + 4(r 



D~<~ C M 



2A 2 



(6) 



This "non-conservation" of currents is originally owing to 
the presence of superconductor condensate which can ab- 
sorb and provide electrons. In the emerged Majorana pic- 
ture, it appears as a consequence of the tunneling-induced 
creation and annihilation of electron pair in the absence 
of bias voltage (across the quantum dot), because of the 
superposition of electron and hole excitations. Based 
on Eq. (^|), one can verify the existence of the Majorana 
fermion from two aspects, (i) Modulating the quantum 
dot level by a gate voltage, the current difference will ex- 
hibit a Lorentzian "lineshape" for the dot level (erj) de- 
pendence, which centers at the Majorana's zero energy. 
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FIG. 2: (Color online) Difference of the steady-state source 
and drain currents, unlike the usual case that is zero, which re- 
veals the Majorana's unique nature of electron and hole (i.e., the 
antiparticle component) excitation. This current difference de- 
pends on the ratio of Tl and Tr, which not only indicates the 
existence (effect) of the Majorana fermion, but also provides an 
access to extract its coupling amplitude to the quantum dot. In 
this plot we set to = 0. 



(ii) As shown in Fig. 2, the current difference depends 
on the asymmetry of Tt and Tr. In particular, we find 
that (only) for = Tt the source and drain currents are 
equal to each other. And, there exists a special setup with 
the right rate Tr equal to TJj, 



transforming the EOM and solving the resultant algebraic 
equations can directly give the result of S a (uj), by noting 
that the integral in Eq. ([)]) is actually a Laplace transfor- 
mation since, sin(wr) = ~(e tWT — e~ 4 " T ). 
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FIG. 3: (Color online) Shot noise which reveals a spectral dip 
and a zero-frequency enhancement effect. The former reflects a 
coherent oscillation in the system which indicates the formation 
of a bound state coupled to the quantum dot, while the latter 
originates from the peculiar nature of the Majorana excitation. 
Ai = A corresponds to coupling to the Majorana bound state, 
while Ai = to a regular bound state. As in Fig. 2, we set also 
eo — in this figure. 
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the current difference reaches a maximum 
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This simple result can also be used to justify the existence 
of the Majorana fermion, and to extract its coupling am- 
plitude (A) to the quantum dot as well. 

Shot Noise. — As recognized in quantum transport, 
the shot noise can carry valuable information beyond 
the steady-state current. We therefore, in the following, 
further explore the shot noise properties of our setup. 
Specifically, we consider the current correlator S a (t) = 
?({5I a (t),6I a (0)}}, with a = L and R and SI a {t) = 
Ia{t) — I'a- Here I* is the steady-state current in the a- 
th lead; and the (quantum statistical) average in defining 
the correlator is over the steady state as well. The shot 
noise spectrum, i.e., the Fourier transform of S a (t), can be 
calculated most conveniently within the number-resolved 
master equation formalism using the McDonald's formula 
11 



S a {w) = 2e 2 uj 



dr i 



(9) 



where (n a (r)) = J2n a n aP (nc "'( T )> and i n l( T )) = 
2 n "■qP'™ q ' ) ( T )- Here, p(" a )(r) is defined by counting 
the electron number "n Q " through the a-th junction start- 
ing from an arbitrary moment in the steady state. There- 
fore, we have (n a (T)) = 7*r. To calculate (n 2 (r)), 
combining the number-resolved master equation, we can 
construct its equation-of-motion (EOM). Then, Laplace 



Figure 3 displays a representative result for the shot 
noise spectrum. First, we notice that a spectral dip ap- 
pears at the frequency us c ~ 2 A (for ed = and small 
cm), which reflects an existence of coherent oscillations 
in the central system, with the characteristic frequency u> c . 
This is an important signature, since it indicates the emer- 
gence of bound states at the ends of the wire with discrete 
energy, gapped from other higher energy continuum. At 
this point, we should keep in mind that the quantum dot is 
coupled to a nanowire, so in the usual case the wire states 
are extended and have almost continuous energies, which 
cannot support coherent oscillations as indicated by the 
spectral dip in Fig. 3. As a comparison, we plot in Fig. 
3 as well the result for the same quantum dot and with 
the same strength (A) coupled to a regular bound state. 
While an "oscillation dip" appears similarly at the same 
frequency, nevertheless the zero (and low) frequency noise 
differs remarkably from the Majorana case. 

To characterize the zero frequency noise, we employ 
the well-known Fano Factor, F a = S a (0)/2el^. For 
symmetric rates = Tl = To, we find Ft = -Ffl and 
denote the Fano factor simply by F. Analytically, we ob- 
tain 
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(10) 



Here F^ R ' is the Fano factor for coupling to a regular (R) 
bound state. Also, in this result we distinguish the cou- 
pling amplitudes A and Ai (introduced previously in the 
coupling Hamiltonian). We notice that A and Ai play iden- 
tical (symmetric) role and the difference, Eq. (10 1, van- 
ishes if any of the amplitudes disappears. We then under- 
stand that the zero-frequency noise enhancement is aris- 
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ing from the peculiar nature of the Majorana excitation. 
Therefore, the noise enhancement effect in Fig. 3 is an- 
other useful signature for Majorana excitation at the ends 
of the nanowire. Also, Eq. (10) can be used to extract the 



important parameters of the Majorana's mutual interac- 
tion (em) and its coupling with the quantum dot (A). In 
particular, based on Eq. (|l0|), an even simpler result can 
be obtained in the limit cm — > by noting Ai = A for 
Majorana fermion. That is, F = 1 - \\V + 2(A/r n ) 2 ]" 1 . 
This result provides a very simple relation between the 
Fano factor and the scaled coupling amplitude (A/To). 

Now we consider Tl i=- Tr. Taking the limit cm — > 
and setting €q = 0, we obtain 
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We notice that the first term in Fl and Fr is the Fano fac- 
tor corresponding to transport through an isolated single- 
level quantum dot [p4||, while the second term arises from 
coupling to Majorana fermion. If Tl ^ Tr, we find 
Fr — Fl = 24(r L -r^)rA 2 /(r 2 + 8A 2 ) 2 ^ 0. (This dif- 
ference vanishes when Tl = Fr.) As a final remark, we 
notice that the steady-state current through the quantum 
dot cannot reveal the Majorana information in the sym- 
metric case (Tl = Fr). However, the Fano factor given 
above carries such information, since the second term in 
F^^R-f does not vanish when y = 1. Note also that, if 
the quantum dot couples to a regular bound state (another 
dot), the zero frequency noise (Fano factor) is the same as 
the first term of the above Fur\, being unaffected by the 
side coupling. 

Experimental Feasibility. — The semiconductor 
nanowire sketched in Fig. 1 can be the InSb wire utilized 



in the recent experiment [|22|). The InSb nanowire has 
a large (/-factor (g ~ 50), and a strong Rashba-type 
spin-orbit interaction (with energy ~50 /ieV). For 
an applied magnetic field of 0.15 Tesla, the Zeeman 
splitting, E z ~ 225 fieV, starts to exceed the induced 
superconducting gap A ~ 200 /icV, thus being able to 
support the emergence of Majorana fermion at the ends 
of the InSb nanowire. Therefore, a low temperature such 
as T = 100 mK can suppress the thermal excitation of 
the Majorana zero mode to higher energy states. As a 
consequence, the condition ksT < eV < E z defines a 
broad range of bias voltage across the quantum dot to 
detect the Majorana's existence and its coherent dynam- 
ics. For the latter purpose, one may tune the tunneling 
rates ^l(r) an d the coupling energy A to the order of 
a few /zeV, which can maintain a coherent oscillation 
between the QD and the Majorana states, provided the 
system has a coherence time longer than nanoseconds. 
Moreover, the above parameters can guarantee a large 
bias condition for electron transport through the quantum 
dot, which is needed to ensure the employed master 
equation applicable. The above estimates indicate that the 
proposed measurement scheme in this work is feasible by 
the present state-of-the-art experiment. 

Summary. — We have proposed a quantum dot trans- 
port scheme to detect the existence and coherent dynam- 
ics of Majorana fermion. To be specific, we considered 
the quantum dot being tunnel-coupled to a semiconductor 
nanowire which is anticipated to generate Majorana exci- 
tations under appropriate conditions. We carried out pe- 
culiar signatures in both the steady state current and shot 
noise for our demonstrating purpose. 
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